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parameters.  Quasisteady  models  are  obtained  for  the  one  and  two 
dimensional  cases,  and  specific  application  is  made  to  both  compres¬ 
sible  and  incompressible  boundary  layer  situations.  Several  analytical 
and  numerical  solutions  are  obtained  for  the  incompressible  models.^ 


Quasisteady  Viscous  Flows 


I  Introduction 


In  an  earlier  paper  [1]  the  author  showed  the  following: 

In  an  internal  region  of  the  flow  where  significant  variations 
occur  only  over  many  mean-free  paths,  the  viscous  terms  of  the 
compressible  Navier-Stokes  Equations  provide  only  a  small  pertur¬ 
bation  on  the  forces  described  by  the  inviscid  equations.  This 
conclusion,  although  consistent  with  known  results  concerning 
strong  shocks  and  thin  boundary  layers ,  presents  difficulties 
in  regard  to  incompressible  flows.  In  [1]  the  following  was 
stated: 

...  there  exist  many  calculations,  with  the  incom¬ 
pressible  Navier-Stokes  Equations,  in  which  the 
viscous  effects  are  pronounced  and  yet  significant 
variations  occur  only  over  many  mean-free  paths. 

The  analysis  of  the  present  paper  indicates  that 
the  use  of  the  time-dependent  incompressible 
Navier-Stokes  equations  in  these  situations  is 
totally  unjustified. 

However,  no  resolution  of  this  problem  was  offered. 

In  later  work  [2,3]  the  author  considered  atmospheric  inviscid 
flows.  The  atmospheric  problem  is  not  unlike  the  incompressible 
boundary  layer  situation:  The  hydrostatic  assumption  in  atmospheric 
models  is  much  like  the  assumption  =  0  of  incompressible  models. 
In  [2,3]  the  approach  was  to  use  quasisteady  analysis  to  obtain 
simplified  mathematical  models  which  are ,  1)  consistent  with  the 
scale  of  the  physical  problem,  and  2)  consistent  with  the  original 
mathematical  model.  The  results  of  these  investigations  leads 
one  to  conclude  that  quasisteady  analysis  is  both  valid  and  useful. 
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In  order  to  use  the  approach  of  the  atmospheric  problem,  it 
is  necessary  that  the  mathematical  model  consist  of  a  system  of 
first  order  hyperbolic  partial  differential  equations.  For  viscous 
flows  such  a  model  exists  in  the  form  of  Cauchy's  Equations  if  an 
appropriate  specification  is  made  for  the  viscous  parameters. 

The  purpose  of  the  present  paper  is  to  apply  quasisteady  analysis 
to  Cauchy's  Equations,  in  which  relaxation  equations  are  used 
to  model  the  viscous  terms. 

Section  II  considers  two  relatively  simple  situations  ,  both 
of  which  will  be  needed  in  the  later  analysis:  1)  a  single  hyper¬ 
bolic  equation  with  a  forcing  term,  and  2)  the  non-Fourier  heat 
conduction  equation.  In  Sections  III  and  IV  quasisteady  models 
are  obtained  from  Cauchy's  Equations  of  motion  for,  respectively, 
the  one- dimensional  and  two-dimensional  cases.  Section  V  considers 
the  problem  of  specifying  boundary  conditions.  In  Section  VI 
application  is  made  to  boundary  layer  situations  for  both  the 
compressible  and  incompressible  cases.  Some  conparison  is  made 
between  these  models  and  existing  models.  Finally,  Section  VII 
discusses  several  analytical  and  numerical  solutions  to  the  in¬ 
compressible  case. 
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Section  II:  Examples  Using  Quasisteady  Analysis 

Two  examples  are  considered  in  this  section.  Both  will  be 
needed  in  the  later  analysis  of  the  viscous  equations. 

A)  A  Single  Hyperbolic  Equation  with  a  forcing  Term 
The  following  equation  is  to  be  considered: 

9t  +  c0x  =  f(0*  *•  t)*  a.1) 


Two  distinct  phenomena  are  involved  in  this  equation: 

(i)  The  hyperbolic  effect  -  through  c,  the  "sound  speed", 
information  is  transmitted  from  the  boundaries  to 
interior  points  of  the  region. 

(ii)  The  forcing  effect  -  the  forcing  term,  f,  does  not 
transmit,  but  rather  "creates"  information  as  a 
point  source. 

Each  effect  occurs  at  its  own  rate.  Our  specific  interest 
was  in  the  consideration  of  quasi-steady  situations  related 
to  these  rates.  One  typical  special  case  of  eq.(l.l)  is  the 
following : 


9-0, 


9t  +  c9x  = 


(1.2) 


where  c,  0q,  and  t-^  are  positive  constants  and  0  <  x  <  L.  The 
characteristic  curves  of  the  hyperbolic  system  have  slope 
^  i  .  Thus,  a  unique  solution  is  determined  by  initial  data 
and  a  boundary  condition  at  x  =  0.  Considering  only  the  boundary 
condition, 


e(t,o>  =  @x(t)  +  0O, 

the  solution  for  t  >  £  can  be  written  as, 


(1.3) 


4 


e(t,x) 


ctr  ei(t  -  !> 


+  9 


(1.4) 


Suppose  first,  that  the  hyperbolic  term  is  dominant.  That 
is,  assume 


ct. 


<<  1. 


(2.1) 


Also,  suppose  the  boundary  data  varies  slowly.  That  is,  assume 
( 0 j_ (t)  -  0^(t  -  £•)  |  is  negligible  for  all  t.  (2.2) 


(This  assumption  was  also  important  in  our  earlier  work, 

[2,  assumption  (2)3).  Then,  from  eqs.  (1.4),  (2.1),  and  (2.2), 

e(t.x)  '  0x(t  -  |)+  e0  ~  e1(t)  +  0Q.  (2.3) 


Returning  to  eq.  (1.2),  the  quasi-steady  solution  is, 


0 


x 


or  0(t,x)  ~  0^(t)  +  0Q. 


(2.4) 


Thus,  assumptions  (2.1)  and  (2.2)  lead  to  a  meaningful  quasi¬ 
steady  solution. 

Suppose  next  that  the  forcing  term  is  dominant.  That  is, 
suppose 

sj”1-  <2-5> 


Then,  from  eq.  (1.4),  0(t,x)  “  0q  except  near  x  =  0;  at  x  =  0 , 
0(t,O)  =  9^(t)  +  9g .  We  see  then  that  the  dominance  of  the 
forcing  term  produces  a  boundary  layer  effect  unless  6^(t)  =  0. 
Furthermore,  in  this  case  a  quasi-steady  solution  would  not  be 
valid  near  the  boundary.  If  0^(t)  =0,  the  solution  is  simply 
©x  =  0,  or  0=  9q.  The  latter  situation  will  arise  in  the  study 


5 


of  the  pressure  equation  in  Section  V. 

B)  The  Non-Fourier  Heat  Conduction  Equation 
The  heat  conduction  equation, 

Tt  =  kTxx*  0  <  x  <  L*  t ->  °*  (3.1) 

involves  the  assumption, 

q  =  -  kTx.  (3.2) 

That  eq.(3.2)  is  a  quasi-steady  equation  is  a  fact  that  was 
heavily  emphasized  by  early  researchers  (see,  for  example, 
Maxwell  [4]).  The  following  non-quasi-steady  equation,  or 
non-Fourier  equation,  has  been  investigated  by  various  people; 
see,  for  example,  refs.  5  and  6; 

Tt  -  -  <x 

(3.3) 

<t  "  -  FJ  <kTx  +  *>• 

In  vector  form  we  obtain: 


Vt  =  AVx  +  BV* 

V  -  (3,  A  =  (_J}„  "b.  B  =  (° 


(3.4) 


-k  ft, 


0  -1/t 


)• 


Eqs.  (3.4)  and  (1.2)  have  the  same  form.  Letting 


c  =  ,  (4.1) 

C1 

the  matrix  A  has  eigenvalues  ±c  with  corresponding  character¬ 
istics  of  slope  ±  i  .  Also,  the  matrix  B  has  the  eigenvalues 

0,  — .  It  follows,  then,  that  both  characteristic  speeds  for 
C1 
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the  hyperbolic  system  have  the  same  magnitude,  namely  i;  however, 
the  forcing  term  always  has  a  zero  eigenvalue,  although  the  other 
may  become  large.  As  in  eq.  (1.2),  it  is  first  necessary  to 
compare  the  magnitudes  of  c  and  t^. 

Eq .  (3.4)  is  very  interesting  and  was  studied  in  some  depth. 
Complete  details  will  not  be  given  here.  A  major  difference 
between  eqs .  (3.4)  and  (1.2)  is  that  in  the  latter  c  and  t^ 
can  be  considered  independent  quantities  while  in  eq.  (3.4)  they 
are  related  through  eq.  (4.1).  There  are  three  time  scales  to 
consider : 


> 

rt 

H* 

II 

rt 

forcing  term  scale, 

(4.2) 

At2  =  Ax/c: 

hyperbolic  scale, 

(4.3) 

o 

At3  =  AX  /k: 

parabolic  scale 

(4.4) 

Ax  is  not  here  intended  as  a  numerical  quantity,  but  rather  as 
some  physically  significant  length.  In  order  to  take  a  quasi¬ 
steady  limit  with  respect  to  the  forcing  term,  we  would  need 

A  t^  =  eAt2,  e  small.  (5.1) 

[Eq.  (5.1)  should  be  compared  with  eq.  (2.5)].  From  eqs.  (4), 

2  2 

we  obtain  At*  =  =  At,  At-.  Then,  from  eq .  (5.1), 

c  k 

At2  =  eAt^  •  (5.2) 

Suppose  a  physical  situation  is  such  that  a  priori  one 
knows  that  At^  is  the  appropriate  time  scale.  If,  then,  e  is 
small,  the  quasi-steady  solution  requires  that  in  eq.  (3.4) 

Vt  =  0,  or  we  obtain  the  steady-state  solution  to  eq.  (3.4). 
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Such  a  solution  can  be  valid;  the  requirement  is  that  boundary 
conditions  for  eq .  (3.1)  vary  slowly  with  time,  in  the  sense 
of  eq.  (2.2) . 

The  general  situation,  however,  would  be  that  e,  of  eq. 
(5.2)  is  small,  but  the  steady-state  solution  is  not  valid. 

This  means  that  it  is  not  possible  to  independently  take  quasi¬ 
steady  limits  of  the  forcing  term  or  the  hyperbolic  term. 
Apparently,  if  e  is  small,  the  forcing  term  effect  occurs  more 
rapidly  than  the  hyperbolic  effect,  but  both  disappear  in  the 
limit  of  e  =  0.  This  situation  is  still  somewhat  unclear.  The 
different  eigenvalues  of  the  matrix  B  are  a  factor  here. 
Analytic  solutions  were  obtained  for  specific  problems,  such  as 
T(t,0)  =  T(t,L)  =  0.  These  solutions  indicate  that  the  limit 
t-^+0  of  the  solution  of  eq.  (3.3)  does  approach  the  solution 
of  eq.  (3.1),  although  convergence  is  somewhat  "precarious". 

The  opposite  case,  e~l,  is  also  surprising,  but  more  infor¬ 
mative.  In  this  case,  since  At^~At2~At^ ,  the  relaxation  effect 
occurs  on  the  same  scale  as  the  other  effects  and  consequently 
cannot  be  deleted. 

The  above  analysis  is  directly  applicable  to  the  Navier- 
Stokes  Equations  in  the  following  ways: 

i)  The  situation  of  e~l  will  correspond  to  a  "mesh  Renolds 
Number"  of  approximately  1. 


ii)  The  two-dimensional  viscous  equations  have  a  "subsystem" 
with  precisely  the  form  of  eq.  (3.4). 
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Section  III:  The  One-Diiaensional  Viscous  Flow  Equations 

The  one-dimensional  equations  of  motion,  Cauchy's  equation, 
[7]  are  as  follows: 


p  +  UP  +  pu  =0, 

t  X  X 


p (u  +  uu  )  +  p  =  a  , 
t  XXX 


(6) 


pc 


(T  +  uT  )  +  pu  =  -q 
v  t  x  x  ^x 


+  Ou 


Where  p  =  pressure,  p  =  density,  u  =  x  component  of  velocity, 

T  =  temperature,  t  =  time  coordinate,  x  =  distance  coordinate, 
cv  =  specific  heat,  q  =  heat  flux,  a  =  shear  stress. 

Upon  introducing  the  quasi-steady  assumptions, 


q  =  -kT  , 
x 


(7.1) 


a  =  pu  , 

x 


(7.2) 


pRT, 


(7.3) 


eq.  (6)  becomes  the  classical  Navier-Stokes  Equations.  In 
the  following,  we  introduce  relaxation  terms  into  eqs .  (7.1) 
and  (7.2): 


p  =  -up  -  pu  , 

t  X  X 


pt  T  pp  1 

u  =  -uu  -  —  T  -  —  p  +  —  o  , 
t  xpxpxpx 


T  =  -uT - —  q  -  u  , 

t  x  pc  ^x  pc  x 

V  V 

q  =  -uq  -  T  -  , 

x  x 


,  U  o 

0  =  — UO  +  —  u  -  — 

t  X  x 


(8) 


p  “  p(P»T). 
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where  k  =  coefficient  of  thermal  conductivity,  y  =  coefficient 
of  shear  viscosity,  and  t^  and  t2  are  relaxation  times. 

In  vector  form, 


Various  features  of  ea.  (9)  were  studied.  First,  the 
eigenvalues  of  A  were  found  to  be  the  solution  of  the 
following  equation  in  X: 
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In  order  that  eq.  (9)  be  hyperbolic  it  is  sufficient 
that  the  roots  of  eq.  (10)  be  real  and  distinct.  A  necessary 
condition,  since  b  >0,  is  that, 

a  >  0.  (10.2) 

A  complete  analysis  of  eq.  (10.2)  would  require  specific 
information  regarding  the  relaxation  times  t^  and  t2-  Although 
such  an  analysis  will  not  be  attempted  here,  eq .  (10.2)  can 
be  considered  a  restriction  on  the  parameters  t^  and  t2 . 
Assuming  eq .  (10.2),  eq.  (9)  will  be  hyperbolic  if, 

a2  -  4b  >  0  (10.3) 


Various  special  cases  were  examined.  For  example,  since 


4b  =  (a  -  2cJ)2  +  4c2c2(1  -  ^ 


(10.4) 


eq.  (10.3)  is  satisfied  in  the  special  cases  where  —  <  1  or 

Cj  =  0 .  Also,  intuitively  there  is  reason  to  suppose  that 

c  ,  c  ,  and  c  are  comparable  quantitities .  If  one  assumes 
s  q  (7 

2  ,  2  a  2,  ,  ,  2  2,,  Ov 

cs  =  Cq  =  Co  =  c  then  a  "  4b  =  (c  '  p  ca)  +  4cqCT(1  '  p>  * 

If  one  further  assumes  that  p  =  pRT  and  that  a  =  0,  then 

a2  -  4b  =  c4[l  +  ^  ]  .  Letting  y  =  14,  the  eigenvalues  are 

found  to  be  X  =  u,  u  ±  1.494c,  u  ±  .8 764c. 

It  is  not  of  course  intended  that  the  above  be  a  defini¬ 


tive  study  of  assumptions  (10.2)  and  (10.3).  Nevertheless, 
at  this  point  it  is  reasonable  to  suppose  that  eq.  (9)  is 
hyperbolic  for  many,  if  not  all,  cases  of  interest.  In 
particular,  this  is  so  for  the  incompressible  limit  c^,  =  0. 

Secondly,  the  normalized  variables  were  considered  [8J. 
This  is  a  somewhat  complex  study,  in  that  the  diagonalizing 
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matrix  for  A  needs  to  be  found  (as  well  as  its  inverse) . 

The  result  is  that  the  normalized  variables  for  eq.  (9) 
are  comparable  to  those  obtained  in  ref.  8  for  the  inviscid 
equations.  This  was  considered  important  since  the  impli¬ 
cation  is  that  perturbation  properties  for  eq .  (9)  are  similar 
to  those  for  the  inviscid  equations. 

Thirdly,  it  is  necessary  to  find  a  set  of  variables  for 
which  the  eigenvalue  corresponding  to  A  =  u  is  in  some  sense 
uncoupled  from  those  involving  the  sound  speed.  Assuming  that 

t„ 


=  constant 


(ID 


a  suitable  set  is, 


Then 


W  =  -A , W  -  B.W, 
t  lx  1 


r*\ 

u  1 

h 

w  = 

T 

Mi 

» 

p*  =  pe  p  a*  =  p  -  o  . 

(12.1) 
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Remark:  If  assumption  (11)  does  not  hold,  it  may  not  be 
possible  in  all  cases  to  achieve  an  uncoupled  form  through  an 
analytic  transformation.  This  does  not  affect  the  validity 
of  the  quasisteady  equations,  but  the  analysis  needs  to  be 
considerably  modified.  It  is  of  interest  to  note  that  the 
term  p*^,  which  in  the  steady-state  limit  reduces  to  pux,  is 
now  a  forcing  term,  and  not  part  of  the  hyperbolic  effect. 

Finally,  the  possibility  of  quasisteady  solutions  was 
considered.  If  a  quasisteady  situation  exists  relative  to 
the  various  "sound"  speeds  and  if  the  flow  velocity  is  small 
relative  to  these  speeds,  then  one  must  use  the  steady-state 
form  of  the  last  four  equations.  The  system  is  then. 


1*  -  p*u  , 

(13.1) 

X  X 

.  -±(p- 

yu  )  , 

(13.2) 

x  p 

x  x 

p-pu 

)u 

-  uT  + 

T 

(13.3) 

pcv  X 

X 

pc  XX 

V 

If  we  then  let  t2-*-0,  e<3-  (13 .1)  becomes  the  continuity  equation. 
Summarizing,  the  momentum  and  energy  equations  are  in  quasi¬ 
steady  equilibrium  with  respect  to  the  continuity  equations. 

In  addition,  a  condition  such  as  eq.  (2.2)  is  also  implied  for 
validity  of  eqs .  (13). 

Suppose  next  that  the  length  of  interest  is  small;  that 
is,  suppose 


Ax  A  =  mean-free  path  =  -  ; 

r  pc 

s 


(U.l) 
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where  Ax  is  again  intended  as  a  physically  significant  length. 
Suppose  also  that  the  relaxation  times  are  on  the  order  of 
the  time  required  to  traverse  a  mean- free  path  at  the  sound 
speed. 


2  _  kCs  k  .  2 

cq  pcvA  (  pcv  )cs' 


(1A.2) 

(14.3) 

(14.4) 


yc 


where,  as  noted  in  [1],  is  essentially  a  Prandtl  number. 

If  this  Prandtl  number  is  near  unity,  then  c  -  c  =  c  =  c, 

q  s  a 

and  the  hyperbolic  time  scale  for  the  problem,  At  ~  is 
comparable  to  the  forcing  time  scale  t^. 

If  u  is  comparable  to  c,  then  in  this  case,  as  in 
the  discussion  of  the  heat  conduction  equation,  eq .  (12) 
should  be  used  as  it  is.  Under  assumptions  (14.1)  and  (14.2), 
eq.  (12)  is  on  the  same  order  of  difficulty,  in  terms  of 
obtaining  a  numerical  solution,  as  the  usual  inviscid  system. 
One  expects  the  solution  to  steady-state  flows  to  be  similar 
to  those  obtained  from  the  Navier-Stokes  Equations.  On  the 
other  hand,  solutions  to  time-dependent  flows  (such  as  shock 
interaction  problems)  could  be  considerably  altered. 

If  | u |  <<  c,  then  eqs .  (13)  could  be  a  feasible  set  of 


quasisteady  equations. 
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Section  IV:  The  Two -Dimens tional  Viscous  Flow  Equations 

The  Cauchy  equations  now  take  the  following  form: 

+  up  +  wp  +  pu  +  pw  =0, 
t  x  z  x  z  ’ 

p (u^  +  UU  +  WU  )  +  P  =  T  +0*. 
t  X  Z  X  z  x’ 

p(w_  +UW  +  ww  )  +  p  =T  +  o**  (15.1) 

txz  zxz  v 

pc  (T  +  uT  +  wT  )  +  p(u  +  w  )  =  -q* 
vt  x  z  xz  x 

-q**+  o*u  +  t (w  +  u  )  +  o**w  , 

Z  X  X  z  z 

where  u,  w  are  velocity  components  in  the  x,z  directions, 
respectively,  q*  and  q**  are  the  components  of  the  heat  flux, 
and  the  2x2  stress  matrix  has  the  form  ,  see  ref.  9, 

p.  46.  Using 


T  =  U(wx  +  uz). 


a*  =  y[2u  -  4(u  +  w  )] 

X  J  X  z 
O**  =  y[2w  -  |(u  +  w  )] 

Z  J  X  z 


(15.2) 


q*  =  -kT  and  q**  =  -kT 
x  z 


Eq.  (15.1)  becomes  the  compressible  Navier-Stokes  Equations. 

As  in  [2]  and  [3],  the  top  of  the  region  will  be  allowed 
to  float  so  as  to  remain  a  streamline.  This  is  accomplished 
by  the  transformation, 
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(16.1) 

(16.2) 

(16.3) 

(16. A) 


where  and  represent,  respectively,  the  position  of  the 
lower  and  upper  boundaries. 

If  we  then  also  add  relaxation  parameters  to  the  viscous 
terms,  eqs .  (15)  become: 

p-  +  up^  +  rp  ^  +  p(u^  +  u^Cx)  +  Pt2w?  =  0,  (17.1) 

p(ut  +  uu^  +  ru?)  +  Pn  +  P?CX  =  AzTj,  +  a*  +  a*^,  (17.2) 

P(wt  +  uwn  +  +  P^z  =  tza**  +  +  tcCx,  (17.3) 

pcv(Tt  +  uTn  +  rV  +  p(un  +  Vx  +  Vz}  = 

-<q*  +  q*Ax)  -  q**Cz  +  a*(Ufi  +  u?cx)  (17. A) 

+  ,<«„  +  w£CK  +  V2>  + 


3E 
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where 


r  "  +  u£x  +  wcz. 


The  variable  propagating  at  the  flow  velocity  is  uncoupled 
from  the  system  by  defining, 


p*  =  pe 


3t„  3t~ 

( — =-  a*  +  — -  a**) 
K  2P  2y  ’ 


0*  =  P  -  o*, 

a**  =  p  -  0** 


(18.1) 

(18.2) 

(18.3) 


With  assumption  (11),  namely  that  —  and  —  are  both  constant, 
eqs .  (17)  become, 


rp*  = 

3p* 

2p 

(a*  +  o**) , 

(19.1) 

+  ru  ) 

=  £  x 

-  0*  -  £  0* 

(19.2) 

£ 

z 

£ 

n  x  £ 

+  rw  ) 

=  T 

+ 

£  T  —  £  0** 

(19.3) 

£ 

n 

x  £  z  £ 

’  )  =  - 

0*U 

0*£  u  -  0**£  w 

£ 

n 

x  £  z  £ 

(19.4) 

;  w  + 

x  £ 

‘zV 

- 

(a*  +  £  q*  +  £  q**) , 
n  x  £  Z  £ 

v  t 


-**  ,  2  /r  2  2.  xp  2 

o?*  +  uo*  +  ro*  =  -pc. u  -  p(£  c. - c„,)u  +  —  c„w 

t  n  £  in  x  l  pT^pTn 


)  rrk 

-  n(c2  -  nr )w? _  n (q* +  ^q* +  ^q**}  +  it 

z  z. 


(19.5) 


TC 

o|*  +  uo**  +  ro**  =  -Pc^u_  -  p(£„c,  - -  4)up  +  -T  c‘w, 


xp  2 

£  " ~3“n  "  KVY3  p  ■  y  ctN 

2  2  (19.6) 


2  ^xTCT  ^CT  n** 

-  pcz<c4  -  nr  )wt "  n  (q* +  cxq*  +  czq**)  +  tt* 

z  3 
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(19.7) 

(19.8) 

(19.9) 

(19.10) 

(19.11) 

Sqs .  (17.5),  (17.8),  and  (17.9)  are  unchanged. 

Suppose  that  solutions  are  desired  in  the  region, 

olr>iL>0l^l1,t>0  (20.1) 

Some  simplification  in  eq.  (19)  is  obtained  by  assuming 
that  the  vertical  height  and  vertical  variation  is  small: 


where 


C1  = 
2 

c2  “ 
2 

c3  = 
2 

c4  " 

2 

Cm 


.  a*  2  4y 

pp  +  T  T  +  3pt7  » 


P 

;**  2 


2 

4y 


Pp  +  —  CT  +  3^t: 


2y 


,  a*  2 

Pp  +  p  CT  Ipt 


.  a**  2  , 

Pp  +  —  CT  + 

PPT 

“5 
p  c 


3 

4y 

3pt, 


|  f2  -  f^|  £  several  hundred  mean-free  paths,  (20.2) 


L  I  (f2 


w 


«  1. 


(20.3) 


From  (20.2)  we  assume  that  vertical  velocities  and  pressure 
gradients  are  small.  In  particular,  we  assume: 


We  therefore  delete  w  £x  from  eq.  (17.5)  and  a*cx 


(20.4) 

(20.5) 

from  eq .  (19 . 2)  . 


gaiggajgjgj;  **-  *  •*****&, 


*0&*sa 


wiwwaiawi 
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The  equations  then  take  the  following  form: 


V  =-A V  -  BV  -  CV, 

t  n  c 


where 


(21) 


V  =  (p*,  w,  T,  q*,  q**,  a**,  a*,  u,  t). 


u 


0 

0 

0 


A  = 


0 

0 


0 

0 

0 


0 

u 

t 

pcv 

0 

0 

2 

tpct 


P 


0 

0 

u 

_k 

C4 

0 

0 

0 

0 

0 


0 

0 

1 

pcv 

u 


0 


P 

0 

0 


0 

0 

0 

0 

u 

0 

0 

0 

0 


0 

0 

0 

0 

0 

u 

0 

0 

0 


p* 

pcv 


pc 


pc. 


V 


B  = 


/' 


0 

0 


\  0 


C  =  ( 


3 (a*  +  a **) 
2  v 


P<\ 


T?x 
C4  X 

k 

s'* 


0 

»ck 


0 

iz_ 

pcv 


r 

2 

pVz 


bl  = 


b3  = 


b5  = 


0*C  -  T  £ 

X  2 

PC 


b2  = 


p(  5xc3 


2  tCzct 


p(  Cxcx 


2  tVt 


p 

pCT^x 

P 

P 

P 

0 

0 

0 

0 

-1  -1 
t  ’  t  ’ 

C4  C5 

zl  i± 
t3*  t2 ' 

o**c  - 
z 

T^x 

pc 

V 

9 

b4  = 

pCZ(c4  ' 

b6  = 

PCz(c2 

p 

0 


^xtct 


0 

0 


), 


c  tcJ 

x  T 


The  solutions  of  I B- XI I  =  0  are  X  =  r,  r,  r  +  /—— 

Ptx  ’ 

,  and  the  solutions  of  |  B^- XI |  =  0,  where  B1  is  the 
indicated  5x5  submatrix.  The  latter  5  are  solutions  of  the 


polynomial 


20 


where 


(r  -  A) [ (r  -  A)  -  a(r  -  A)  +  b] 


(22) 


a  =  —  b.  + 


c  + 


p  4  pcvt  ^  z  pcvt4  x 


4..  . 


b  =  (  — ~  r  +  -—V-  4  )(b.  ~  -  -----  ) 

pcvt4  z  f>V5  x  4  p  P 


The  discussion  of  the  one-dimensional  system  given  in  Section  III 
is  applicable  to  eq.  (22).  In  particular,  we  assume  the  following: 


i)  The  solutions  of  the  4th  degree  polynomial  in  eq.(22) 
are 


A  =  r  ±  2^,  r  ±  c 2> 


ii)  | r |  <<  c. 


where  e,  =  min(e, ,  c„,  ) . 

h  Li  pt^ 


(23.1) 

(23.2) 


iii)  As  a  consequence  of  assumption  (20.2),  the  relaxation 

times  are  comparable  to  the  hyperbolic  time  scale,  ^3  3) 

(f2  -  f1)/ch. 


In  the  discussion  of  boundary  conditions  we  will  need  the 
correspondence  between  eigenvalues  and  the  variables.  These  are 
as  follows: 


i) 

A  = 

r,  r,  r  correspond  to  p*. 

0*,  and  q*, 

(24.1) 

ii) 

A  = 

AT 

r  ±  y— —  correspond  to  u 
ptl 

and  t  , 

(24.2) 

iii) 

A  = 

r  ±  2^,  r  ±  22  correspond 

to  w,  T,  q**,  and  a**. 

(24.3) 

Remark:  Eqs.  (24)  are  not  quite  accurate  in  that  the  variables 

involved  in  the  above  matrix  need  to  be  modified  so  as  to 
uncouple  the  variable  corresponding  to  A  =  r.  However,  because 
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of  assumption  (20.3),  £  is  small  and  the  appropriate  variables 
will  be  small  perturbations  of  the  present  ones.  The  result 
will  not  affect  the  final  quasisteady  equations  and  so  the 
derivation  is  not  carried  out  here. 

The  eigenvalues  of  the  matrix  A  of  eq.  (21)  are  given 
by, 


Al|  =  (u-A)3[(u-A)2 


Pti 


][(u-A)^  -  a(u-A)2  +  b]. 


(25) 


where 


2  * 
ci 


v  4 


b  = 


pc  t 

v  4 


a*  c2 

,2  T  , 

(c  i - —  ). 


The  solutions  of  the  4th  degree  polynomial  in  eq.  (25) 


are  easily  seen  to  be  real.  Consequently,  the  eigenvalues  of 


A  are  (A  =  u,u,u,  u  ±  / 


r~ 


pti 


u  ±  c^,  u  ±  where  c^  and 


c^  correspond  to  the  4th  degree  polynomial.  It  is  assumed 
that , 


<  =  min (c 


(26) 


The  correspondence  between  eigenvalues  and  variables  is  as 
follows : 


ij 

ii) 

iii) 


A  =  u,  u,  u  correspond  to  p* ,  q**,  a**. 


A  =  u  ±  correspond  to  w  and  t, 


pt 


A  =  u  ±  u  ±  correspond  to  T,  q*,  0*,  u. 


(27.1) 

(27.2) 


(27.3) 
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As  in  the  atmospheric  problem,  the  goal  here  is  to  con¬ 
sider  quasisteady  situations  induced  by  "distortion"  in  the 
spatial  scale;  that  is,  we  consider  boundary  layer  situations. 
The  physical  scale  assumptions  for  the  present  problem  can 
be  stated  as  follows  [2,  assumption  2] 

i)  The  physical  region,  given  by  eq.  (20.1),  ,9R  , 
is  such  that  C f 2  _  f^)/^  <<  ^/c^  ■  v  •  ; 

ii)  Boundary  conditions  and  initial  conditions 
are  such  that  the  flow  variables  will 
experience  significant  variations  only  (28.2) 
over  time  scales  which  are  large  compared 
to  (f2  -  f1)/ch. 

Remark:  Assumptions  (28)  would  be  typically  satisfied  in 
thin  boundary  layers.  In  many  such  problems  the  vertical 
mesh  has  caused  severe  numerical  problems  [10,11]. 

Because  of  assumption  (28.1),  we  need  not  consider  the 
possibility  of  simultaneously  having  quasisteady  limits  with 
respect  to  both  q  and  C.  By  assumption  (28.2),  however, 
quasisteady  limits  with  respect  to  £  are  feasible.  We 
conclude,  therefore,  as  follows:  By  eqs .  (24.2)  and  (24.3) 
and  assumption  (23.2),  the  variables  [u , t ,w,T , q** , a**}  are 
in  quasisteady  equilibrium  with  respect  to  p* ;  by  assumption 
(23.3),  the  variables  q*,o*  are  in  quasisteady  equilibrium 
with  respect  to  p*.  The  viscous  parameters  are  then, 


T  =  p(wn  +  t  Ur), 

n  z  c 


°*  =Iu(un  +  W  "!  uCzV 


a**  =  p(|w^z  -f  un  -f  W’ 


(29.1) 

(29.2) 


(29.3) 


23 


where  in  eq.  (30.1)  it  has  been  assumed  that  p*-*p. 

Summarizing,  the  quasisteady  equations  (30)  are  proposed 
as  an  approximation  to  eqs.  (17)  if  the  several  assumptions 
(20.2,  20.3,  20.4,  20.5,  23.1,  23.2,  23.3,  28.1,  28.2}  are 
satisfied. 
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Section  V:  Boundary  Equations  for  the  Quasisteady  Model 

The  approach  used  to  obtain  boundary  equations  for  the 
quasisteady  model  will  be  the  same  as  that  described  in  [3]: 
First,  at  each  boundary  an  appropriate  set  of  characteristic 
equations  and  boundary  conditions  is  specified  for  the  time- 
dependent  model;  secondly,  we  attempt  to  simultaneously  take 
the  quasisteady  limit  of  the  internal  partial  differential 
equations  and  the  boundary  equations. 

As  a  specific  problem,  we  will  consider  viscous  flow 
over  a  flat  plate:  n  =  0  will  represent  a  position  somewhere 
downstream  of  the  leading  edge  and  r,  =  0  represents  the  wall. 

The  nine  equations  needed  at  each  boundary  for  the  time- 
dependent  model,  eq.  (21),  will  consist  of  boundary  conditions 
(whose  number  is  determined  by  the  number  of  incoming  character 
istics)  plus  characteristic  (linearized)  equations. 

We  assume  that, 

n  =  0  is  an  inflow  boundary  and  q  =  1  is 
an  outflow  boundary  (Except  at  x,  =  0 ,  u  0  (31) 

at  n  =  0  and  n  =  1)  . 

From  eqs .  (25),  (26)  and  30,  at  n  =  0  six  eigenvalues  are 
positive,  indicating  six  incoming  characteristics,  and  three 
are  negative.  The  nine  equations  can  be  chosen  as 


p  =  h1(t,  O  , 
q**  =  h2(t,rj  , 
o**  =  h3(t,rj  , 
w  =  h4(t,r,)  , 
T  =  h5(t,  O  , 
u  =  h  e  ( t ,  C  )  , 


(32.1) 

\ 

1  (32.2) 

J 
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where  h^  are  specified  functions,  plus  the  three  character¬ 
istic  equations  which  will  have  the  form  (see  [3]  for  more 
detail) : 


allwt  +  a12Tt  =  a13 

a21Tt  +  a22q£  +  a23CT|  +  a24ut  =  a25 

a31Tt  +  a32qE  +  a335g  +  a34ut  =  a35 


(32.3) 


At  n  =  1  there  are  three  incoming  characteristics.  The 
equations  can  be  chosen  as  follows: 


w  =  0 , 

(33.1) 

q*  =  0, 

(33.2) 

Pt  -  “ut. 

(33.3) 

(see  [3]  for  a  further  discussion  of  (33.3)  as  an  outflow 

lateral  boundary  condition.  The  six  characteristic  equations 

will  consist  of  the  partial  differential  equations  for  p*, 

t 

q^*  Og*  and  three  equations  of  the  form  (32.3). 

At  c,  =  0,  r  =  0  or,  referring  to  eqs .  (22)  and  (23),  the 
following  equations  can  be  chosen:  p*,  a*,  and  q*  can  be 
calculated  from  their  differential  equations;  the  boundary 
conditions  can  be  taken  as, 


u 


w 


fi  /ii  +  <c)2'’ 


(34.1) 


T  =  h7(t,n)  , 


5S: 
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and  the  characteristic  equations  have  the  form, 


a41u£  +  a42x£  a43 ’ 

a51Tt  +  a52wf  +  a53  *¥  +  a545!*  =  a55-  (34-2> 

a61Tt  +  a62wt  +  a63’’E*  +  a64°l*  "  a65  • 


At  ?  =  1,  we  again  have  r  =  0.  Thus,  the  equations  for  o* , 
a*,  and  q*  can  be  used.  The  characteristic  equations  will 
have  the  form  (34.2)  and  the  boundary  conditions  are  chosen 
as  , 


x  =  0, 


q**  =  0  =5>T^  =  0. 


Pt 


ux 


(34.3) 


(34.4) 


(  See  C3J  for  a  discussion  of  the  last  equation  in  the 
context  of  an  upper  boundary  condition.) 

Eqs .  (24)  were  obtained  by  assuming  that  the  variables 
{u,  w,  T,  x,  q* ,  q**,  a*,  a**}  were  in  quasisteady  equilib¬ 
rium  with  respect  to  p* .  Quasisteady  boundary  conditions  are 
obtained  by  applying  the  following  principles  [2,  eqs.  (15 
and  (16)  ]: 


If  a  boundary  equation  involves  any  one  of  {u~, 

Wt’  ^t’  Tt’  and  if  the  equation 

is  driven  by  flow  conditions  internal  to  the  region 
of  computation,  then  the  equation  will  be  put  in 
quasisteady  equilibrium. 


] 

i 


1 


(35.1) 
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The  basic  quasisteady  assumptions,  (28.1)  and 

(28.2)  require  that  eqs .  (29 . 1) - (29 . 4)  and  eqs. 

(30 . 2)  - (30 . 4)  be  applied  at  all  points  of  the  (35.2) 

lateral  boundaries,  except  possibly  at  the  corner 
points . 

The  following  is  thereby  obtained: 

a)  At  n  =  0,  eqs.  (32.2)  are  deleted  because  of  (35.2) 
and  eqs.  (32.3)  are  deleted  by  (35.1).  Thus,  p  is 
obtained  from  eq.  (32.1)  and  the  remaining  variables 
are  determined  by  quasisteady  equations. 

b)  At  n  =  1 ,  eqs.  (33)  and  all  characteristic  equations, 
except  the  equation  for  p*,  are  deleted  because  of 

(35.2).  Thus,  p  is  obtained  from  eq.  (30.1)  and 

the  remaining  variables  are  determined  by  quasi¬ 
steady  equations . 

c)  At  the  upper  and  lower  boundaries,  C  =  0  and  1, 
the  equations  for  o,v  and  q*,  as  well  as  eqs.  (34.2) 
are  deleted  because  of  (35.1).  Thus,  p  is  obtained 
from  eq .  (30.1)  at  both  C  =  0  and  1.  Also  eqs.  (34.1) 
are  retained  at  x,  =  0  and  eqs.  (34.3)  are  retained 

at  c  =  1.  (the  last  equation  in  (34.3)  is  retained 
because  it  is  intended  as  a  simulation  of  external 
flow  conditions  [3]). 

Summarizing,  the  proposed  quasisteady  model  is  as  follows: 

i)  p  is  obtained  at  n  =  0  from  eq.  (32.1)  and  at 
all  other  points  from  eq.  (30.1). 

ii)  u  is  obtained  from  eq.  (30.2)  and  the  boundary 
conditions  (34.1)  at  c  =  0  and  t  =  0  at  s  =  1 . 

iii)  w  is  obtained  from  eq.  (30.3),  the  boundary  con¬ 
ditions  (34.1)  at  t,  =  0,  and  a  condition,  to  be  discussed 
below,  at  x,  -  1. 
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iv)  T  is  obtained  form  eq.  (30. A)  and  the  boundary 
conditions  T  =  h(t,n),  where  h  is  specified,  at  c,  =  0  and 
T  =  0  at  e  =  1. 

v)  Eq .  (16.4)  is  a  first  order  hyperbolic  partial 
differential  equation  in  the  variable  f2  (f^  is  specified) 
and,  because  of  eq.  (30),  requires  a  boundary  condition  at 
n  =  0.  We  will  impose, 

(f2)£  =  0  at  n =  0.  (36) 

The  third  condition  of  eqs .  (34.3),  p  =  u  ,  has  not 

cs 

been  used.  It  is  the  missing  condition  for  w.  Its  use 
requires  an  additional  iteration  between  eqs.  (16.4),  (34.4), 
and  the  calculation  of  u  given  above. 

Remark:  Clearly,  eqs.  (30,2),  (30.3),  and  (30.4)  are  coupled 

and  would  in  fact  need  to  be  solved  simultaneously.  The  above 
description  of  the  model  is  intended  to  display  the  primary 
dependence  between  the  variables  and  the  differential  equations. 

Finally,  we  need  to  consider  the  corner  points.  As 

discussed  briefly  in  [3],  inconsistencies  can  arise  here  in 

regard  to  boundary  conditions: 

i)  4=0,  n  =  0:  p  is  specified  by  eq .  (32.1),  T 

is  specified  or  h^(t,0)  =  hy(t,0)  -  see  eqs.  (32.2)  and  (34.1), 

u  and  w  are  specified  by  eqs.  (34.1). 

ii)  c  =  1,  n  =  0:  p  is  specified  by  eq.  (32.1)  and 

YPut 

the  remaining  conditions  are  T  =  0,  t  =  0,  and  pt  =  —  2~- 

^  s 

iii)  f,  =  0,  z,  =  1:  p  is  calculated  from  eq.  (30.1),  and 
the  other  variables  are  as  at  c  =  0,  n  =  0. 

iv)  c  =  1,  n  =  1:  pis  calculated  from  eq.  (30.1),  and 
the  other  variables  are  as  at  c  *  1,  n  =  0. 


Section  VI:  Boundary  Layer  Approximations 

Assumption  (28.1)  has  the  effect  of  requiring  that  the 
width  of  the  region  be  small  relative  to  the  length.  The 
following  assumption  insures  that  the  length  itself  is  large: 

L  is  large  relative  to  the  mean- free  path,  so  that 
significant  variations  in  the  x-direction  occur  only  (37) 
over  many  mean-free  paths. 

From  the  results  in  [1]  we  conclude  that  the  viscous  terms, 
in  the  n-direction ,  of  the  Navier-Stokes  Equations  can  be 
deleted.  Assuming  this  conclusion  also  holds  for  the  Cauchy 
Equations,  eqs .  (30)  reduce  to  the  following: 


pt 

+  upn  + 

rpC 

_ 

-  p  (u  +  r  u, 

n  x  ^ 

+  4 

zV  - 

(38. 

.1) 

0 

-  -p<uun 

+ 

ruc 

>  +  4*,ur.c  * 

V 

(38. 

■  2) 

0 

=  -p(uw 

n 

+ 

rw^ 

)  + 

-  c 

zpc 

(38. 

.3) 

+ 

^zp 

(  3  r-z”r,  ■  3 

cxu 

cv 

0 

=  -pcv(uT 

+ 

n 

rT 

c>  -  P<u„  + 

r,  u 

x  c 

+ 

+ 

2  2  4  2 

(38. 

4) 

PC 

X  +  7  V 

where,  on  the  basis  of  assumptions  (20.3)  and  (20.4),  several 
viscous  terms  involving  c,  have  been  deleted  from  eq.  (38.4). 
Remark:  Assumption  (37)  could  have  been  applied  before  making 
the  transformation  to  {t,  n,  c).  In  this  case  the  viscous 
terms  involving  cx  would  not  appear  in  eq.  (38.3).  This  is 
consistent  with  the  well-known  result  in  boundary  layer 
theory  that  the  final  approximation  is  dependent  to  a  certain 
extent  on  the  coordinate  system. 
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Eqs.  (38)  are  similar  to  Stokes'  Equations  for  compressible 
flow  [12,  eqs.  (20)],  which  are  obtained  as  the  limit  of  the 
Reynolds  Humber  approaching  zero.  Letting, 


M  =  Mach  Number  =  —  , 

cs 

(39.1) 

:e  =  Reynolds  Number  =  , 

(39.2) 

A  =  mean-free  path  =  —  , 

^  /•» 

(39.3) 

c 


we  obtain, 

Re=M(^).  (39.4) 

L/A  is  large,  by  assumption  (37),  but  there  is  as  yet  no 
restriction  on  M  (assumption  (26)  was  introduced  only  for  the 
purpose  of  specifying  boundary  conditions).  Thus,  eqs.  (38) 
are  valid  for  large  Reynolds  Number,  as  long  as  the  scale 
assumptions  are  satisfied. 

Eqs.  (38)  should  also  be  compared  with  the  boundary 
equations  commonly  used  for  time  dependent  compressible  flows, 
as,  for  example,  given  by  Lagerstrom  [12,  eqs.  (20~12)  ]  .  Very 
fundamental  differences  are  seen: 

i)  Lagerstrom' s  equations  retain  the  time  derivative 
in  the  horizontal  momentum  equation  and  the  energy  equation, 
whereas  both  have  been  deleted  in  eq.  (38). 

ii)  Lagerstrom' s  equations  have  replaced  the  vertical 
momentum  equation  with  the  condition  p^  =  0 .  In  eqs.  (38) 
a  further  assumption,  to  be  discussed  below,  is  required 
before  such  a  simplification  can  be  attempted. 

An  incompressible  version  of  eqs.  (38)  is  obtained  by 


assuming  P^  =  0 .  Eq.  (38.4)  is  then  uncoupled  from  the 
remaining  equations.  With  the  notation, 


p  =  p(p)  and  p  =  c  , 


(40.1) 


eqs.  (38)  become, 


pt  +  upn  +  rp?  =  “pc  (un  +  cxuc  +  ?z wc>’  C40 . 2) 

p(uun  +  ru^)  +  pn  =  0,  (40.3) 

p(uV  4zu(3W  iV <40'4> 


Eqs.  (40)  can  be  compared  with,  for  example,  Prandtl's 
boundary  layer  equations  [9,  eqs.  (7.7)  &  (7.8)].  The  essential 
differences  would  be  as  follows: 

i)  In  Prandtl's  equations,  eqs.  (40.2)  and  (40.4)  will 
be  replaced  respectively  by  u^  +  txu^  +  =  0  and  p^.  =  0. 

ii)  In  Prandtl's  equations,  the  time  derivative  term 
would  be  added  to  eq .  (40.3). 

Boundary  equations  for  eqs.  (40),  corresponding  to 
those  given  in  Section  V,  would  be  as  follovs: 


4  =  0: 

u- v£- 

Cx  2 

w  =  f,  /[l  +  (yV],  eq .  (40 . 2)  ; 
t  ^z 

(41.1) 

4=1: 

ur  -  0,  Pt  = 

^  u  ,  eq.  (40.2) ; 
s  c 

(41.2) 

ii 

o 

p  specified, 

eqs.  (40.3)  and  (40.4), 

(41.3) 

S3 

II 

eqs .  (40 . 1) , 

(40.2),  (40.3). 

(41.4) 

.  A  further  simplification  of  eqs.  (40)  can  be  achieved 
by  analyzing  the  time  scales  associated  with  these  equations. 
Since,  in  the  time  dependent  form,  the  right  side  of  eq.(40.1) 
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is  a  forcing  term,  then  with  our  assumptions  the  hyperbolic 
time  scale  for  eq.(40.1)  is, 


Atl  u 


(42.1) 


Integrating  eq.(40.4)  in  the  form. 


^  -  -  IT  uc  +  taxV  +  +  bl- 


where 


3  [r> 

bl  "  zr^-  o  p< 


uw^  +  rw  )d^, 


(42.2) 


p*  =  p(t,  n,  1) , 


and  substituting  into  eq.(40.4),  we  obtain, 


PE  +  PPn  +  ^  +  b2' 


where 


b2  =~pc  (un  +  +  V' 

Ah  - 

t2'l^  ' 


Next,  defining  w£  and  w*  by, 


w  =  w*  +  w 


and 


r  w*  = 
4zwC 


-<un  +  cxu«>- 


eq .  (40.2)  becomes, 


p£  +  upn  +  rp ^  =  -pc  £zwc 


(42.3) 


(43.1) 


(43.2) 


(43.3) 


Integrating  eq.  (43.3)  approximately,  with  the  time  increment 
At-p  and  then  differentiating  with  respect  to  c,  we  obtain, 


(Pn  “  (P  )c  -  At1[(upn  +  rp^)  ]c-  At1C2(pc  w^)c. 

where  the  superscripts  n  and  n+1  denote  respectively  times 
t  and  t  +  At^ .  With  eq.(43.4),  eq.(40.4)  becomes, 


/  9  A 1 1 

Wz'1  + 


4  2 


■  3'cxczyucc 


+  f^z^xVc  +  p(uwn  +  -  Ati?z(pc2)cw^ 

+  Cz[pn  -  Atx(upn  +  rp ^ ) n ] £ . 

At^  3  L  1 

From  eqs.(39),  =  jf(j)  (^)  .  Assuming  that  the  Mach  Number 

is  not  large,  say, 


M  <  1, 


we  conclude,  from  assumption  (37),  that 


>>  1. 


It  is  therefore  appropriate  (see  Section  II. A)  to  take  the 
quasisteady  limit  of  eq.(42.3),  which  thereby  becomes, 


p  =  p* ,  or  p  s  0 . 


At  the  same  time  eq.(43.5)  becomes, 


4  2  *  1. 

inr  w*  - 


2  2  2 
+  ^MCZ(CX)CU;  +  P(uw„  +  rwr )  -  At1?,(pc  )rw 


l^z'^  '  £  e 


“  Atl«zucV 


Because  we  have  taken  an  additional  quasisteady  limit, 


boundary  eqs.(41)  need  to  be  reexamined  in  terms  of  (35.1) 
and  (35.2).  Eq.(40.2),  at  4  =  0  and  t,  -  1,  presents  some 

difficulty:  It  is  assumed  that  the  equation  reflects  only 

internal  properties  of  the  fluid  at  c  =  0,  while  at  z,  =  1  i 
gives  information  regarding  the  external  flow  (since  it  is 
applied  at  the  dividing  streamline  between  the  boundary 
layer  and  the  external  flow).  Thus,  eq .  (40.2)  is  deleted 
at  c  =  0  but  retained  at  c  =  1.  The  proposed  model  is  then 
as  follows: 

i)  At  {0  <  r,  <  1,  0  <  n  <  1), 

eqs . (40 . 3) ,  (43.2),  (45.1),  (45.2). 

ii)  At  { c  =  0 ,  0<  n<  1), 


w*  =  f,  /[l  +  (%2]  , 
t  ^z 

w  =  0 . 

e 

iii)  At  { C  =  1,  0  <  n  £  1}, 


u^  =  0, 

pt '  ?v 


w. 


pc 


2  (?t  +  up„> 


n  ’ 


f?  =  w  +  u  — -  , 
t 


iv)  At  k  “  1,  n  =  0  }, 


(46.1) 


(46.2) 


(46.3) 


(46.4) 


p*  specified, 
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pt  ■  uf 


w. 


pc 


(Pt  +  upn)  • 


Remark:  The  possibility  of  inconsistencies  arising  at  corner 

points  was  discussed  in  [3].  In  eqs .  (46)  such  a  difficulty 
occurs  at  {^  =  1,  n  =  0}:  Since  p*  is  specified,  either 
u  =  0  or  p_  =  ^  u.  must  be  deleted.  This  could  result  in 
some  type  of  "inlet  adjustment"  length  in  the  solution. 

Finally,  it  is  clear,  from  eq .  (44.2)  that  the  w*^  term 
in  eq.  (45.2)  is  negligible.  In  addition,  under  usual  con¬ 
ditions  in  which  incompressible  models  are  applied,  one  expects 
the  remaining  terms  on  the  right  side  of  eq.(45.2)  to  be  small. 
In  any  case,  with  the  assumption, 


|w£ |  <<  w*,  (47) 

w£  becomes  a  purely  diagnostic  quantity.  The  resulting  model 
in  (p,  u,  w}  is  as  follows: 

Eqs.  (46)  with  w  =  w*  and  w£  deleted.  (48) 

In  comparing  model  (48)  with  Prandtl's  boundary  layer 
equations  [9,  eqs. (7.7)  and  (7.8)],  one  sees  that  the  remaining 
discrepancy  is  in  the  time  dependence.  In  eqs.  (48),  time 
dependence  occurs  only  through  boundary  conditions,  while  in 
Prandtl's  equations,  time  dependence  occurs  through  the  time 
derivative  term  in  the  horizontal  momentum  equation.  In  some 
problems,  particularly  those  involving  small  perturbations  in 
time,  this  difference  might  be  significant. 


Summarizing,  we  have  derived  several  models  for  possible 


use  in  boundary  layer  problems: 

i)  Eqs .  (38)  are  proposed  for  compressible  flow. 

However,  this  model  is  still  incomplete: 

a)  The  derivation  needs  to  be  carried  out  for  the 

t  O 

case  of  variable  — (assumption  (11)), 

b)  The  effect  of  the  forcing  terms  needs  to  be 
considered  in  eq.(38.1). 

ii)  Eqs. (46)  and  (48)  are  proposed  for  incompressible 

boundary  layers.  Solutions  to  these  equations  will 
be  presented  in  the  next  section. 

Although  the  above  models  are  intended  to  be  time- 
dependent,  a  serious  question  remains  concerning  the  types  of 
time-dependent  behavior  retained  in  quasisteady  models .  The 
following  remarks  address  this  question. 

Remark  1:  Quasisteady  models  assume  that  appropriate  time 
dependence  for  the  quasisteady  variables  is  provided  by  boundary 
conditions.  A  clear  example,  discussed  in  [2],  is  the  following 
Suppose  the  flow  is  completely  one-dimensional  (say,  in 
the  n-direction) ,  but  a  thin  c~region  is  maintained.  Then, 
eqs.  (38)  might  still  be  valid,  but  clearly  cannot  by  themselves 
provide  time-dependent  propagation  of  the  flow.  The  important 
point  here  is  that  eqs.  (38)  are  internal  equations,  while  the 
boundary  conditions  are  considered  separately.  Appropriate 
boundary  conditions  could  be  to  apply  the  full  one-dimenstional 
equations  at,  say,  q  =  0 ;  these  equations  would  be  unaffected 
by  the  internal  quasisteady  procedure  with  respect  to  c, . 

In  this  regard  note  that  eqs. (46)  and  (48)  were  derived  for 
specific  boundary  conditions. 
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Remark  2:  As  a  second  example,  consider  flow  in  a  narrow 
channel  bounded  by  two  parallel  flat  plates.  Without  going 
through  all  details,  the  quasisteady  model  corresponding  to 
eqs .  (48)  would  be  as  follows: 

i)  At  the  "internal"  points  {0  <  c  <  1,  0  <  n  <  D. 
the  same  set  of  equations  would  be  used:  eqs.  (40.3),  (43.2), 
(45.1),  (45.2), 

ii)  At  { c  =  0  and  r,  =  1,  0  <  n  £  U,  u  =  0,  w  =  0, 

2 

P £  =  -ec  V 

iii)  p*  specified  at  the  inlet  line  n  =  0 . 

(Note  that  in  this  problem  f~  =0  and  the  equation 

t 

p  =  -^  ut  is  not  appropriate) . 

It  is  not  expected  that  the  above  model  can  predict  ♦'he  develop¬ 
ment  of  the  flow,  beginning  with  stationary  flow  at  time 
zero  and  a  time-dependent  pressure  input  at  n  =  0.  This  can 
be  explained  by  noting  that  two  of  the  major  assumptions  are 
not  satisfied  for  this  problem: 

a)  One  expects  the  initial  propagation  for  this 
problem,  at  any  point  in  the  flow,  to  resemble  a  one-dimensional 
compression  wave  (following  which,  the  vertical  viscous  terms 
dominate  the  flow).  Thus,  assumption  (37)  is  not  valid  at 

all  times. 

b)  Unless  the  region  is  extremely  thin  (say,  20 
mean-free  paths),  assumption  (28.1)  will  not  be  satisfied 
in  the  neighborhood  of  the  moving  compression  wave  (here 
we  interpret  L  to  be  the  width  of  the  wave) . 

Remark  3:  Eqs.  (46)  cannot  predict  the  formation  of  the 
boundary  layer  resulting  from  impulsive  motion;  assumption 


(28.2)  would  not  be  valid.  They  should,  however,  be  applicable 
to  some  types  of  perturbations  of  an  existing  boundary  layer. 

An  example  of  such  a  problem  will  be  given  in  Section  VII. 
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Section  VII:  Solutions  to  the  Incompressible  Boundary 
Layer  Model 

In  Part  A  we  consider  analytic  results  regarding  eqs .  (46), 
and  in  Part  B  we  present  some  numerical  solutions  to  several  time- 
dependent  and  steady-state  problems. 

A)  Analytic  Considerations 

Assuming , 

i)  that  w^  ~  0,  (49.1) 

ii)  that  the  convective  terms  in  eq.(40.3)  are  (49.2) 

negligible , 

iii)  that  -^=pc  can  be  considered  constant  in  eq .  (49.3) 

c 


(46.3) , 

eqs.  (46)  take  the  following  form: 

=  f2pr),u(t,  n,0)  =  u^(t,n,l)  =  0,  (50.1) 

w  =  -fu  +  (fi  +  4f  )u  ,  w(t,n,0)  =  0,  (50.2) 

^  n  n  ^ 

pc  =  p0c0utoP,  pe(t,0)  specified,  (50.3) 

f2  =  wtop  -  f2  utop,  fo(t,0)  ==  f 2 (0 , n)  =  fQ,  (50.4) 

zt  n  ^  “ 

Mi.  II 

where  f  =  f2  -  f^,  Pe  =  p(t,n,l)  -  pQ ,  P  denotes  4=1. 

It  has  been  assumed  that  u  =  0  at  t  =  0  and  that  f,  =  0. 

t 

From  eqs . (50 . 1)  and  (50.2), 

f2 

u  =  -  (ut0p)(42  -  20,  ut0p  =  -  ,  (51.1) 


r 
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w  =  -  c2)f(utop)  -(t2-2c)f1  utoP-(^! 


From  eqs.(50.3)  and  (51.1), 


ut°P  “  2y  P0C0(U  P)n* 


Q2)fr)UtOTj.  (51.2) 


or 


(utoP)  =  -  ^  ut0P. 
n  f  2 


From  eqs .  (50.4)  and  (51.2), 

= 


§  (fut0p)„. 


n 


From  eqs. (51.3)  and  (51.4), 
fo  = 


0  top 

3  -  2i>' 


(51.3) 


(51.4) 


(51.5) 


From  eqs. (51.3)  and  (50.3),  one  finds  the  initial  conditions 
to  be  as  follows: 


utop  (0  ,  n)  =  e 


=  „toP  p"  2Arl/f 


0 


0. 


u 


top  1 


p0co  e 


Pc (t , 0)  . 


(52) 


The  steady-state  solution  must  satisfy 

(53.1) 

(53.2) 

Note  that  eq.(53.1)  implies  the  conservation  of  mass  condition, 
(j^1fudi)n=  0.  From  eqs.  (53.2),  (51.3),  and  (50.3)  the  steady- 
state  solution  is. 


(futop)  =  0 

n 

and 

ff  =  2A . 
n 


41 


(The  last  inequality  assumed  that  f,  =  0)  . 

in 

It  would  appear  from  eq.(52)  that,  by  specifying  pe(0,0)  =  0, 
the  solution  could  start  with  stationary  flow  and  generate  a 
boundary  layer.  However,  because  fg  >  0  must  be  arbitrarily 
specified,  this  procedure  would  not  seem  to  be  reasonable. 

Eq.  (53.3)  indicates  that,  as  expected,  the  steady-state 
boundary  layer  grows  as  /x.  Also,  eq.(51.5)  indicates  that  the 
rate  of  growth  of  the  time-dependent  boundary  layer  would  be  ^-utoP 
(assuming  X  <  f ) . 

Note  in  the  above  that  the  model  is  "consistent"  in  terms 
of  mass  conservation  independently  of  eq.(50.3).  Thus,  as  is 
commonly  done  in  boundary  layer  problems,  one  could  specify  p^ 
arbitrarily  and  obtain  solutions.  However,  because  c  =  1  is 
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always  a  streamline,  not  all  functions  for  p  will  give  sensible 

n 

answers.  For  example,  p  =  constant  will  give  uto^  -  constant 
and,  from  eq.(52.1),  f  =  constant. 

The  above  indicates  that  defining  the  edge  of  the  boundary 
layer  as  a  streamline  is  considerably  different  than  assuming 
that  the  "outer  flow"  is  achieved  at  the  edge  (see,  for  example, 

[9,  eqs. (7.5) , (7.7) , (7.8) ]) .  The  use  of  the  streamline  requires, 
for  a  complete  interpretation,  the  existence  of  an  additional 
transition  region  from  the  edge  of  the  boundary  layer  to  the  free 
stream  conditions.  On  the  other  hand,  the  streamline  definition 
avoids  the  inconsistency  in  other  boundary  layer  models  of 
nonzero  vertical  velocity  at  the  edge  of  the  boundary  layer. 

Inconsistencies  at  the  outer  edge  were  also  discussed  by 
Telionis  and  Gupta  [13],  where  it  was  stated  that  u,  p,  and  T 
should  satisfy  the  one-dimenstional  time-dependent  inviscid 
equations  at  the  outer  edge.  Thus,  only  one  of  the  three 
variables  is  arbitrary.  Likewise,  in  the  quasisteady  model  one 
upper  boundary  condition  is  specified,  but  the  other  variables 
are  calculated  from  quasisteady  equations  rather  than  the  time- 
dependent  inviscid  equations. 

B)  Numerical  Solutions  to  Eqs.  (46) 

The  problem  is  that  of  flow  over  a  flat  plate.  The  initial 
stationary  fluid  was  chosen  with  properties, 

p0  =  3.2632  x  10'4  ^  (54.1) 

in 


Tq  =  260  °K . 


(54.2) 
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(These  are  typical  atmospheric  values  at  a  height  of  60  km. [14]). 

We  then  have , 

c  =  c0  =  /yRT0  =  323 . 24m/ sec ,  (54.3) 

u  =  2.556  X  10_6/rTT5"  =  4.121  x  10"5  ,  (54.4) 

X  =  ~r~  =  -3907  x  10-3  m,  (54.5) 

P0C0 

P  =  P0  -  P0RTq  =  24.35  -k-S  2  ,  (54.5) 

m  sec 

As  discussed  in  Section  VI,  there  is  an  initialization 
problem  at  time  zero.  Our  procedure  was  as  follows: 

i)  An  initial  f9(n,0)  =  f0  and  an  initial  pressure 

Z  Zq 

p(0,0)  =  Pg  +  p*  are  specified. 

ii)  With  the  exception  that  the  quantities  f9  and  f,  are 

zt  t,- 

maintained  at  zero,  the  remaining  equations  of  eqs.(46)  are 
integrated  timewise  until  a  steady-state  is  achieved.  These 
values  are  then  the  initial  conditions  for  the  complete  time- 
dependent  problem,  and,  for  small  Mach  number,  should  be  in 
close  agreement  with  eqs.(52). 

In  the  first  problem  to  be  discussed  we  chose, 


f9  =  .02m  ~  51 . 2 A , 

(55.1) 

L  =  2.5m  ~  6399 \ , 

(55.2) 

i — < 

o 

ii 

o 

(X 

1 

o 

o 
' — ' 

cu 

(55.3) 

f1(t,n)  =  0. 

(55.4) 

The  above  lengths  are  consistent  with  the  scale  assumptions  of 
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the  model.  The  pressure  value  is  small  so  that  the  computation 
will  represent  a  low  Mach  number. 

One  expects  the  time  step  in  the  problem  to  be  determined 
by  the  time-dependent  partial  differential  equations.  Upon 
examining  eqs .  (46),  one  sees  that  only  one  such  equation 
remains,  namely  that  for  f-2: 

f2  =  ■  uf2  +w.  (56. 1) 

t  n 

Eq.(56.1)  is  a  simple  first  order  hyperbolic  partial  differential 
equation.  The  time  step  restriction,  the  stability  condition 
for  the  upstream  differencing  scheme,  is. 

At  <  4J1  .  (56.2) 

The  first  computation  was  chosen  with  a  5  x  13  spatial  mesh, 

A;  =  .25  and  An  =  L/12,  or 


AC  =  .25 

An  =  .2083  (57) 

At  =  2. 


(As  will  be  seen  in  the  results,  eq.(55.3)  produces  a  maximum 
u  of  .09482,  which  then  accounts  for  the  choice  of  At).  The 
computations  were  then  repeated  with  the  mesh  cut  by  a  factor  of 
2.  All  results  discussed  below  are  for  the  finer  mesh. 

Figures  1-6  show  results  obtained  for  this  problem.  Figure 
1  shows  the  initial  p  and  u  distribution  at  the  top  after  the 
initialization  has  been  completed.  The  pressure  drops  from 
.01  at  n  =  0  to  .000066  at  n  =  1,  while  u  drops  from  .09482 
to  .000627. 
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Figure  2  shows  the  time-dependent  development  of  the  boundary 
layer,  assuming  the  above  initial  conditions.  Figure  3  displays 
the  boundary  layer  growth  at  two  specific  points.  Although  the 
precise  rate  of  growth  of  the  boundary  layer  is  difficult  to 
recover  from  the  numerical  results.  Figure  4  attempts  to  show 
the  rate  of  growth:  At  each  point  x,  we  plot  the  time  at  which 
the  boundary  layer  achieves  907.  of  its  final  steady-state 
thickness.  Figure  5  shows  the  development  of  the  pressure 
field  and  Figure  6  shows  the  final  steady-state  horizontal 
velocity  distribution. 

The  initial  and  final  steady-state  values  are  in  good 
agreement  with  the  results  established  in  Part  A,  as  can  be  seen 
from  the  following  tabulation. 


utoP(0,2. 

5) 

top 

Pc 

(0,2 

.5) 

utop(0,  1 

•  25) 

top 

Pc 

(0,1 

.25) 

f  (°°>  2 

•  5) 

f(°°,l 

.25) 

utoP< 

°°>  2 . 

5) 

uC°P( 

°°,1. 

25) 

wtoP( 

00 , 2 . 

5) 

wtoP< 

°°,1- 

25) 

f  (°°,2.! 

j)utop(“,2 

Calculated 


.6270  x 

10'3 

.6613  x 

10“4 

.7711  x 

10"2 

.8133  x 

10'3 

.06817 

.04958 

.02773 

.03792 

.3659  x 

10“3 

.6624  x 

10"3 

Analytic 


.7176  x 

10~3 

.7569  x 

h- * 

o 

1 

-c> 

.8249  x 

10-2 

.8700  x 

10~3 

.06562 

.04851 

.02889 

.03838 

.3440  x 

10'3 

.6132  x 

10"3 

► 

} 

► 

J 


eq .  (52) 


eq. (53.3) 


eq. (53.4) 

eq. (53.6) 


.001890 


Table  I 


.001897 


eq. (53.1) 
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The  results  are  generally  as  expected.  The  following  com¬ 
ments  are  perhaps  worth  noting: 

i)  w  is  not  displayed,  but  throughout  |w  |  <  10”^|w|. 

k-  c 

Thus,  eq.(43.2)  is  the  appropriate  equation  for  vertical  velocity. 

-  2 

ii)  In  this  calculation,  |wj  ~  10  |u|.  Nevertheless, 

because  of  eq.(46.3),  the  vertical  velocity  is  a  significant 
quantity . 

iii)  As  noted  earlier,  it  is  not  claimed  that  the  above 
results  represent  the  correct  time  development  of  a  real  boundary 
layer  (the  difficulty  being  in  the  initialization).  However,  the 
time  scale,  particularly  as  exhibited  in  Figure  4,  seems  reason¬ 
able:  The  time  rate  of  growth  seems  to  be  on  the  order  of  the 

horizontal  velocity. 

iv)  The  input  quantity  fg  is  clearly  related  to  the  distance 
one  must  travel  from  the  leading  edge  before  the  quasisteady 
assumptions  are  valid. 

The  accuracy  of  the  computations  was  checked  as  follows. 

Since  the  numerical  scheme  is  first  order  in  the  time-dependent 
equations,  the  error  in  the  computation  should  be  reduced  by  at 
least  a  factor  of  two  if  the  mesh  is  reduced  by  a  factor  of  two. 

Let  V(t,  n ,  4 ,  h)  represent  the  quantity  V  obtained  at  (t ,  n ,  4) 

with  a  mesh  {At,  At,  An),  and  V*(t,  q,  t)  represent  the  correct 
value.  Letting  then  j  represent  the  mesh  ^-1,  we  can 

make  the  following  two  checks  [2]: 

i)  A2  =  ! V(t ,  n,  4,  £)  -  V*(t,  n,  4) |  <  j  Ax  = 

|V(t,  n,  4,  h)  -  V*(t,  n,  4)  |  , 


□ 


1 

I 


(58.1) 
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ii)'  A2  =  |V(t,  n,  4,  £)  -  V(t,  n,  4,  j)  |  £  \  Ax  = 


i  | V(t ,  n,  4,  &)  -  V(t,  n,  4.  h) j . 


(58. 


With  h  representing  the  mesh  given  by  eq .  (57)  and  with  V* 
values  taken  from  the  analytic  column  of  Table  1,  we  can  tabulate 
as  follows : 


V 

A2 

A1 

a2/ax 

utop(0,2.5) 

.906  x  10"4 

1.958  x  10"4 

.463 

f(“,  2.5) 

.00255 

.00598 

.426 

f  (°°.  1  •  25) 

.00087 

.00243 

.358 

utop(cc,2.5) 

.00116 

.00242 

.479 

utop  (oo>  i .  25) 

.00046 

.00188 

.245 

wtop(°o,2 . 5) 

.00002190 

.00004708 

.465 

Table  II 

One  concludes  that  the  test  defined  by  eq.(58.1)  is  satisfied. 

To  test  eq.(58.2),  a  third  calculation,  using  a  mesh  defined 
by  h/4,  was  made.  The  results  tabulated  below,  again  indicate 
that  this  test  is  also  satisfied. 


V 

A2 

A1 

A2/A 

f(100. 

2.5) 

.00012 

.00044 

.273 

f(100. 

1.25) 

.00054 

.00148 

.365 

u(100 , 

2.5) 

.00072 

.00145 

.49  7 

u(100. 

1.25) 

.00064 

.00136 

.471 

w(100 , 

2.5) 

.46  x  10'5 

1.63  x 

10-5 

.282 

w(100 , 

1.25) 

1.76  x  10-5 

6.98  x 

’O'5 

.252 

p(100, 

2.5) 

.76  x  10'4 

1.55  x 

1 

o 

r-f 

.490 

p(ioo, 

1.25) 

.73  x  10'4 

1.41  x 

10'4 

.518 

A 


Table  III 


.518 
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Having  achieved  a  steady-state  solution,  the  model  should  be 
able  to  predict  the  time-dependent  effect  of  subsequent  pertur¬ 
bations  of  the  flow  field.  In  the  next  problem,  the  lower  boundary 
(namely,  the  flat  plate)  is  perturbed  to  an  "angle  of  attack" 
position  by  the  following  equation: 


1 6  n  f 


max 


=  < 


(trto) 


(t  -  t0)2(t  -  tp2 


t  <  tg  and  t  > 


tg  <  t  C 


tg  +  ti 


t0+tl 


(59) 


Eq.(59)  is  such  that  f^  is  linear  in  n  and  achieves  its 
maximum  f^  =  f-j™3*  at  n  =  1  and  t  =  — y— ' .  For  the  sample 
calculation  we  chose, 

fimax  =  25  =  (,i)L(  (60.3 

tg  =  130  sec.  and  t-^  =  170  sec.  (60. 

With  this  input  the  plate  will  rotate  from  zero  incidence  at 
t  =  130  sec.  to  5.71°  at  t  =  150  sec. 

Results  of  the  above  computation  are  shown  in  figures  7  and 
8.  At  time  130  sec.,  the  flow  field  is  near  steady-state.  The 
plate  then  rotates  at  an  average  velocity  of  .0125  m/sec.  at 
n  =  1  during  the  next  20  sec.  Thus  the  vertical  velocity  is 
comparable  to  the  horizontal  velocity  during  a  portion  of  this 
period.  Figure  7  shows  that  the  boundary  layer  thickness  first 
expands  and  then  experiences  a  slow  decay.  The  final  steady-state 
value,  f  =  .0667,  occurs  near  t  =  270  sec.,  and  represents 
approximately  a  2 . 87»  decrease  in  the  boundary  layer  thickness. 
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Figure  8  shows  horizontal  velocity.  One  sees  an  initial 
drop  followed  by  a  rapid  overshoot,  and  then  a  gradual  decay. 
Steady-state  values  are  achieved  more  rapidly  then  the  thickness 
shown  in  figure  7.  The  final  value,  u  =  .02820  represents  a 
2.470  increase. 

The  steady-state  vertical  velocity  for  this  problem  is  on 
the  order  of  107o  of  the  horizontal  velocity,  while  we  still  have 
lwel  ~  10'4|w/'I|  . 

Remark:  The  actual  "length  of  computation"  along  the  physical 


2 


2.5125m. 
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Figure  1:  Initial  u  and  p  distribution 


Figure  2:  Boundary  Layer  Growth 


Figure  3:  Boundary  Layer  Growth  at  Two  P 


Figure  4: 


Figure  5:  Pressure  Distribution 


Figure  6:  Horizontal 
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Boundary  Layer  Thickness  at  Nonzero  Angle  of  Incidence 


Figure  8:  Horizontal  Velocity  after  dlate  Rotation  at 
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